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Abstract 



We investigate the combinatorial complexity of geodesic Voronoi diagrams on poly- 
hedral terrains using a probabilistic analysis. Aronov et al. [AdBTOS] prove that, if one 
makes certain realistic input assumptions on the terrain, this complexity is Q{n-\-m^/n) 
in the worst case, where n denotes the number of triangles that define the terrain and 
Q ' m denotes the number of Voronoi sites. We prove that under a relaxed set of assump- 

tions the Voronoi diagram has expected complexity 0{n + m), given that the sites have 
a uniform distribution on the domain of the terrain (or the surface of the terrain). 
^ ' Furthermore, we present a worst-case construction of a terrain which implies a lower 

■ bound of r2(nm^/'^) on the expected worst-case complexity if these assumptions on the 
\^ I terrain are dropped. As an additional result, we can show that the expected fatness of 

a cell in a random planar Voronoi diagram is bounded by a constant. 

(N 

1 Introduction 

Voronoi diagrams on terrains are a basic geometric data-structure that have a variety of 
^ I applications in many areas: geographic information science (gis) jAJTOll IDG081 IPA061 

■ ICGlOj . robot motion planning jTS89j . mesh generation |KWR97j and image analysis |SF04t 
lWDB+08j to name a few. The geodesic Voronoi diagram of point sites on a polyhedral 
terrain is a subdivision of the surface into cells, according to the set of sites, such that every 
cell contains exactly the surface points which are closest to the site that is associated with 
the cell. Here, the distance is measured by the length of the shortest path on the terrain. It 
is tempting to believe that in practice - that is, given that the terrain is well-behaved - the 
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complexity of such a geodesic Voronoi diagram should be linear, because of its similarity to 
the Euclidean Voronoi diagram of point sites in the plane. 

However, in the worst case, this complexity can be much higher, even if one makes certain 
realistic assumptions on the shape of the terrain. Indeed, Aronov et al. |AdBT08] show that 
the worst-case complexity is 0(n + m^fn) for a certain class of well-behaved terrains, where 
n is the number triangles that define the terrain and m is the number of Voronoi sites. This 
shows that assuming realistic input indeed brings the complexity down, but it is still far 
from being linear. They conjecture that, in order to prove a linear bound, one needs to make 
further assumptions on how the sites are distributed. Our purpose in this paper is to study 
the complexity of a geodesic Voronoi diagram if we assume that the sites are being chosen 
randomly from the terrain. 



Previous work. Analyzing the expected complexity of geometric structures for random 
inputs has a long history in computational geometry. See for instance the work of Renyi 
and Sulanke jRS63] and Raynaud |Ray70| on the complexity of convex hulls of random 
points. Weil and Wieacker give an overview of related results in jWW93] . Naturally, Voronoi 
diagrams (and their counterpart, Delaunay triangulations) have been analyzed in this way 
as they are a fundamental data-structure, used in fields such as mesh generation |Rup95] , 
surface reconstruction [Deyll] , molecular biology, and many others. It is well-known that the 
size of the Voronoi diagram of point sites in IR^ is quadratic in the worst-case, however it is 
near-linear in most practical situations. To address this dichotomy, people have investigated 
the complexity when the point sites are: (i) generated by a random processes, (ii) well spaced, 
(iii) have bounded spread, or (iv) were sampled from surfaces according to curvature. See 
|DEG08J and references therein for more information on such work. In particular, there 
is a vast amount of work on Poisson Voronoi Diagrams (PVD). Here, the domain has a 
density associated with it (say, the area). The probability of n points to appear in an 
area of measure /i has, as the name suggests, a Poisson distribution parameterized by the 
area. Similarly, the distribution of points selected into disjoint areas is independent. Poisson 
Voronoi diagrams are used in many areas, such as physics, biology, animal ecology, and 
others. See |OBSC00| IJN04] and references therein. However, this work does not seem to 
have considered geodesies at all. 

In this paper, we are interested in the complexity of geodesic Voronoi diagrams on poly- 
hedral terrains. Moet et al. |MvKvdS08j were the first to study this complexity using a set 
of parameterized assumptions that describe realistic terrains. In this approach, one assumes 
that a certain property, for example, the maximum slope of the terrain, can be bounded by 
a constant independent of the input size. This allows one to avoid certain worst-case config- 
urations which are highly unlikely to occur in practice. Instead, the analysis is confined to 
classes of well-behaved inputs and consequently this method is described as using realistic 
input models. Moet also did an experimental validation of the used parameters |Moe08] 
and confirmed that the parameters indeed behave like constants on realistic terrains. The 
realistic input models introduced in this work have also been adopted by subsequent papers. 
As such, Aronov et al. |AdBT08] improved the bounds given by Moet et al. and showed 
that (i) the bisector between two sites has worst-case complexity 0(n), (where n denotes the 
number of triangles of the terrain) if the triangulation is low density and the lifted triangles 
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have bounded slope; and (ii) that the worst-case combinatorial complexity of the Voronoi 
diagram is Q{n + rriy/n), (where m denotes the number of sites) if in addition the triangles 
are of similar size and the aspect ratio of the domain is bounded. The realistic assumptions 
made in these papers are described in more detail in the next section. 

Finally, note that Schreiber and Sharir [SS08] showed how to compute an implicit rep- 
resentation of the geodesic Voronoi diagram on the surface of a convex polyhedron, in time 
0((n + m) log (n + m)), so that the site closest to a query point can be reported in time 
0(log (n + m)). Schreiber [Sch07j also extended their method for single-source shortest paths 
to the case of non-convex polyhedra using several realistic input models. Naturally, these 
analyses do not inform about the complexity of the explicit Voronoi diagram. 

Our results. We study the expected complexity of geodesic Voronoi diagrams on terrains. 
To this end, we use realistic input assumptions on the terrain, and sample the sites uniformly 
at random from the domain of the terrain. See Section[2]for the exact definitions. In Section[3] 
we show that under these assumptions the complexity of the geodesic Voronoi diagram is 
indeed linear. That is, we show that the complexity is bounded by 0{n + m), where n is 
the complexity of the terrain, and m is the number of sites being randomly picked. The 
constants in the asymptotic analysis depend on how well-behaved the terrain is, which is 
formalized using the input models described in the next section. See Theorem 13.31 for the 
exact result. In Section H] we analyze the expected complexity if these assumptions on the 
shape of the terrain are dropped. In particular, in Theorem 14.121 we show a lower bound of 
fli^nm?/^). This lower bound, in a sense, justifies the input assumptions made previously, 
since it implies that the randomness assumption by itself is not sufficient if we want the 
geodesic Voronoi diagram to have a low complexity. The construction that leads to this 
lower bound is intricate and requires a careful balancing of the variance of the distances of 
the sampled sites, and how closely they can be packed together. Furthermore, in Appendix 1X1 
we show that in expectation the Voronoi cells generated by a uniform sample in the plane 
are in expectation fat; that is, they are nicely behaved in some sense. 

Organization. In Section [2] we introduce the concepts used. Section [3] contains the poofs 
of the upper bound on the complexity of the Voronoi diagram. The lower bound is proved 
in Section |H Finally, we conclude in Section [51 The expected fatness of a Voronoi cell is 
analyzed in Appendix 1X1 

2 Preliminaries 

2.1 Voronoi Diagrams on Terrains 

A polyhedral terrain T is defined by a triangulation A of n vertices V in IR^, a convex 
domain D C IR^ which contains V, and a height function on these vertices. The surface of 
the terrain is defined by the triangles of A lifted according to this height function. We refer 
to T simply as a terrain and we denote the set of edges of the triangulation with E. For 
simplicity of exposition we restrict our discussion to the case where D is the unit square, 
however, our results can be easily extended to the more general case of convex regions with 
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bounded aspect ratio. For two points q, s G D, we denote their euclidean distance in the 
{x, j/)-plane with ||q — s||. When q and s are hfted to the surface of T, we define their geodesic 
distance to be the length of the shortest path connecting them that is constrained to he in 
the surface of T, and we denote this value by da-(q, s). 

The (geodesic) Voronoi diagram of a set of m points on T (which are called sites) 
is a planar subdivision of the surface of 7, where every cell of the subdivision is associated 
with exactly one site, and such that for any point in the cell the associated site is the closest 
site, where the distances are measured using the geodesic distance. We denote the Voronoi 
diagram with Vor(P), where P denotes the set of sites, and we call a cell of the subdivision 
a Voronoi cell. The bisector between two sites q and s on the surface of 7 is defined as 
the set of points p, such that p has the same distance to q and s. The Voronoi diagram can 
be represented as the structured set of curves and straight lines which delineate the Voronoi 
cells and which are subsets of the bisectors between these points. We call a point which is 
incident to at least three cells a Voronoi vertex and we call each maximally connected 
subset of the bisector incident to two Voronoi cells a Voronoi edge (note that two cells 
can have multiple edges between them). Usually one assumes general position of the sites 
so that no two sites are equidistant from a terrain vertex, which ensures that bisectors are 
one-dimensional and that the Voronoi cells subdivide the terrain surface without overlap, 
see also |AdBT08] . In our case the sites are randomly sampled, and so general position is 
implied. 

Since a terrain T is defined by a height function over a domain D, there is a natural 
bijection between points of 7 and points of D. Hence, the various objects defined in the 
previous paragraph can be viewed either in T or in D. Generally in the paper we shall refer 
to these objects by their projection in D, unless otherwise stated. 

2.2 Input Model 

The main idea of realistic input models is to parametrize certain properties of the input, 
which are suspected to capture contrived configurations leading to high complexities or 
running times. In cases where there exists a high discrepancy between the theoretical bounds 
and the complexities observed in practice, it is often useful to analyze the complexities not 
only as a function of the input size, but also with respect to these parameters. This sometimes 
leads to more informative asymptotic bounds. As such, the realistic input assumptions do not 
only distinguish between "good" and "bad" input, instead they enable a more differentiated 
view on which inputs are 'better' or 'worse'. 

In this paper, we use the following realistic input model. A set of line segments is X-low 
density if and only if the number of edges that intersect an arbitrary ball, which are longer 
than the radius of the ball, is smaller than A. Low density has been used in the analysis of 
many different geometric problems, see |dBKSV02] for an overview. 

To model a realistic terrain we adopt the realistic assumptions made in |AdBT08] . 
According to these assumptions, there exist constants A and ^ independent of n, such that 

(i) the set of edges of the triangulation, is a \-low density set, and 

(ii) any line segment embedded in the lifted triangulation has slope at most ^. 

We now state some useful facts that follow from these assumptions. For notational ease 
in the rest of the paper we define the constant /3 = ^/T+~^. 
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First, the number of pairs of objects from two low density sets, which intersect each other 
is hnear in the total number of objects in these sets. This can be easily verified, the idea is 
to charge each intersecting pair to the smaller of the two objects. 



Fact 2.1 ( |AdBT08] ) Let A he a set of n objects with X-low density and let B be a set of 



m objects with (p-low density, then the number of pairs of objects {u,v) E A x B, such that 
u intersects v is 0{Xm + (pn). 

Second, since the slope is bounded by a constant, the geodesic distance is the same as 
the euclidean distance up to a constant factor, and similarly, geodesic disks have an area 
that is approximately the same as that of planar disks in this case. 

Fact 2.2 ( [AdBTOS] ) For any two points q,s E D, we have that ||q — s|| < d'j(q,s) < 

/3||q-s||. 

Lemma 2.3 Let D be a geodesic disk of radius r on the surface of a terrain with bounded 
slope C, and let A denote its area. We have that 7r(r//3)^ ^ A < nf^r"^. 

Proof: Let c be the center of D and let Cp be its projection. Let -Dr//3 and Dr be the planar 
disks with center Cp and radius r/(3 and r, respectively. Let T^/^ and denote the portion 
of the terrain that lies directly above -D,,/^ and D^, respectively. 

Clearly the projection of D is contained in Dr and so if we can bound the area of Tj. 
then this will bound the area of D. Observe that consists of a set of triangles from A, 
which have been clipped at the boundary of Dr- It is a well-known fact that if we lift any 
such (clipped) triangle up to the terrain then its area can increase by at most a factor of /3. 
Therefore the total of Tr is at most 7r/3r . 

Now consider the disk -Dr-//?- First observe that Tr/^ must have area at least as large as 
that of Dr/ff. Second, note that D must contain T,,/^ (since by Fact l2.2l the distance between 
any point T^/^ and c is at most (5{r / (5) = r). Therefore the area of D is at least 7r(r//3)^. ■ 



2.3 Complexity of the Voronoi Diagram 

The complexity of the Voronoi diagram is measured by the complexity of the structured set 
of curves and line segments that delineate the Voronoi cells. This set consists of pieces of 
bisectors and it can be characterized as follows. Again, we adopt the definitions used in 
[AdBTOS] . 

For most of the points on a bisector, the shortest path to either site will be unique. If the 
shortest path is not unique, we call p a breakpoint. The breakpoints partition the bisector 
into a set of curved pieces which we call chords. The combinatorial complexity of the 
Voronoi diagram is now defined as the sum of (i) the number of Voronoi vertices, (ii) the 
number of breakpoints of Voronoi edges, and (iii) the number of intersections of the chords 
of Voronoi edges with the triangulation of the terrain. 

We continue with some useful facts and lemmas used in the analysis of the complexity. 
First, it was observed by Moet et al. that the number of breakpoints of the Voronoi diagram 
is bounded by n, since each of them can be attributed to a terrain vertex. To see why this is 
true, imagine walking along the bisector while sweeping the shortest paths to either site from 
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the current position. The points on the bisector, where the edge sequence of the shortest 
path changes, are the breakpoints. 

Fact 2.4 ( |MvKvdS08] ) Given a terrain 7 which is defined by a triangulation with n 
vertices, the number of breakpoints of any Voronoi diagram on T is smaller than or equal to 
n. 

Furthermore, we will use the following result by Aronov et al. [AdBTOS] . 

Lemma 2.5 ( [AdBTOS] ) Given two points q and s, the set of chords that form the bisector 
of q and s on 7 (projected to the {x,y) -plane) is 0{^)-low density. 

We remark that Aronov et al. use this result to show that the bisector has linear com- 
plexity. However, note that this result does not imply that the overall set of chords of the 
Voronoi diagram is low density, which would imply a linear complexity for the whole dia- 
gram. Consider for example the situation, where all the sites lie close to each other on a 
straight line, and all the triangles of the terrain surface are coplanar. In this example, the 
bisectors are pairwise parallel lines, which extend from one side of the domain to the other 
and could therefore lead to a quadratic complexity Voronoi diagram by intersecting many 
triangles of the terrain. 

Finally, we observe that the number of Voronoi vertices and edges is linear in the number 
of sites, as the following lemma and corollary testify. This fact was also observed by Aronov 
et al.., we include an independent proof for the sake of completeness. 

Lemma 2.6 Let 7 be a terrain and let P be a set of m points. Then the number of Voronoi 
edges and Voronoi vertices ofVor{P) is 0{m). 

Proof: For m < 2 the claim is clearly true, since we have at most one bisector, which 
contributes exactly one Voronoi edge. For m > 2, we argue as follows. 

First, observe that the cells in this Voronoi diagram are connected. Indeed, consider a 
point p that belongs to the interior of the cell of s G P. Consider the shortest path vr from p 
to s, and consider any point q G vr. If q is closer to some other site t than to s, then we have 
that 

dT(p, s) = dT(p, q) + dT(q, s) > da(p, q) + da(q, t) > dT(p, t) , 

but this is a contradiction to p being in the interior of the cell of s. 

Now, consider the dual graph 9 of the graph formed by the Voronoi vertices and Voronoi 
edges. In this graph, every vertex corresponds to a Voronoi cell and every face corresponds 
to a Voronoi vertex. Note that we can derive a geometric embedding of this graph by using 
the sites as vertices and picking an arbitrary point on each Voronoi edge and connecting it 
by its shortest path to either site to form a (possibly curved) edge between two vertices. 

It is well-known that a cell in the Voronoi diagram might not be simply connected. 
Indeed, consider a mountain surrounded by a plane. If we place a site s on the top of the 
mountain, and a site t at the bottom of the mountain (and the mountain slope is large 
enough) then the Voronoi cell of s would be completely surrounded by the cell of t and thus 
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would create a 'hole' in this cell. In S, the two vertices that correspond to the cells of s and 
t would be connected by an edge, which is incident to only one face in S. 

Furthermore, it is known that the dual graph can have multiple edges between two sites. 
To see this, again, place s on the top of the mountain and place two sites t and r at the 
bottom, such that the bisector between t and r intersects the mountain. The boundary 
between the cells of t and r would contain two Voronoi edges from the same bisector. 

However, the dual graph S is planar and connected and as such its Euler characteristic 
is 2. Hence v — e + f = 2, where e denotes the number of edges, / the number of faces and 
V the number of vertices of S- 

Now, by definition, every Voronoi vertex is incident to at least three Voronoi cells. There- 
fore, every face of S is incident to at least three edges of S. Since S is planar, every edge of 
9 is incident to at most two faces of S- Hence, we have that 3/ < 2e, which implies that 
3/ < 2(w + / - 2), and therefore it holds that / < 2u - 2 = 2m - 2. It follows that the 
number of Voronoi vertices is in 0{m). Applying Euler 's formula again, we obtain that also 
the number of Voronoi edges is in 0{m). m 

Corollary 2.7 Let 7 be a terrain and let P be a set ofm points. Let n be a sub-square which 
intersects k Voronoi cells in their projection. The number of Voronoi edges which intersect 
Vor(P) n □ m their projection is in 0{k). 



3 Upper bound 

We prove the following lemma first in the planar case and then extend it to terrains with 
bounded slope. The bounded expected complexity then follows by examining the number of 
intersections of the chords with the terrain triangulation in an i/m x ^/m grid. 

Lemma 3.1 Let P be a set of m points, sampled uniformly at random from a unit square, 
and let n be a sub-square contained in the unit square of side length \j ^fm. Then the expected 
number of points in P that contribute to 'Vor(P) fl □ zs 0(1). 

Proof: We place a sequence of exponentially growing disks centered at the center point of 
□. Let ri = ^2\ for i = 0, . . . , A; = [ig (i.e. ro is the radius of the circumscribed 

circle of □). Let dj be the disk of radius r^, which is clipped to the unit square and let 
% = di\ di_i, for z = 1, . . . , /c. 

Let the points in P be labeled pi, . . . , p^. Observe that the expected 
number of points from P that fall into d2 is (Trrj)^ = = Svr = 

0(1). Hence we do not need to worry about their contribution to 
Vor(P) n □. Otherwise, we claim that a point pj which falls into "Jli 
for i > 2, can only contribute to Vor{P) fl □ if dj_2 contains no points 
of P. Assume for the sake of contradiction that the Voronoi cell of Pj 
intersects □, and there exists some point p^ in P that lies in dj_2. By 
construction, we know pj has distance greater than (rj_i —Tq) to any point in do and p^j has 
distance at most (rj_2 + ^^o) to any point in do. Hence we have that d(pj, do) > Tj^i — = 
2ri_2 — tq > ri_2 + ro > d(p„, di) for i > 2 and as such every point in do is strictly closer to 
Pu than Pj, where d(p, X) = miuqgx ||p — q||- 
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Hence it is sufficient to bound the expected number of points Pj which fall into an annulus 
% such that d,_2 is empty. For pj we define the indicator variable X- which is equal to 1 
if and only if pj G 3?,, and the indicator variable which is equal to 1 if and only if no 
other point falls into dj_2. Hence a given point pj can contribute to Vor(P) fl □ if and only 
if XjY/ = 1 for some value of i. Now, we know that 

Pr\Xi = 1] < nr^ - = ^(2^ - 2^(-^)) = ^4-\ 

and 

Pr[Y,^ = 1] < (^1 - ^^rl,j < exp[--rl,{m - l)j = exp(^ A_ Lj < exp(-4^-=^) 

where a factor of 1/4 was added in the bound for Y/ due to boundary effects that might 
arise from the position of □ in the unit square. Hence the number of points that can affect 
'Vor(P) n □ is bounded by J2j J2i>2-^i'^i ^ which in expectation we have, 



E 



j i>2 j i>2 i>2 



by linearity of expectation and the independence of and Y/ for all i and j. ■ 

Lemma 3.2 Let 7 be a terrain with bounded slope C,- Let P be a random sample ofm points, 
either sampled uniformly from T, or uniformly from the unit square and then lifted vertically 
up to 7. Let n be a sub-square contained in the unit square of side length 1 / ^/rn. Then the 
expected number of points in P that contribute to the portion of the Voronoi diagram of P on 
7 that lies above □ is 0{f3^). 

Proof: Follows by a careful adaptation of the proof of Lemma [3.11 To accommodate for the 
larger distances on the surface, we increase the radii of the disks slightly. 

So, let Ti = -^={2(3 + 1)*, for i = 1, . . . , fc. Let d, be the disk (in the plane) of radius rj 
centered at the center of □. And let = dj \ dj_i, for z = 1, . . . , A;, as defined above. For 
points p, q G dj_2 and s e = dj \ dj_i, we have that 

dT(p,q) < /3||P - q|| < /32ri_2 < r-i^i - ri_2 < ||q - s|| < da(q,s) . 

As such, as before, for a point in P fl 3lj to affect the Voronoi diagram on □ requires that 
dj_2 is empty of any points of P. 

Now, consider the case that the points are sampled from the terrain. Define Xl and 
as in Lemma 13. 1[ Note that the area of the terrain can only increase from one by lifting the 
individual triangles. We have that 

^ v^^i 1 area on terrain of 3ii ^/ 9 o \ 7r/3 ^ ,n,- , ^ xq,- 9\ 

Pr X/ = 1 < : < aijxrl - Trr^i < -^ 2/3 + if ^ - 2/3 + if ^"^ 

area of terram 2m 

<^(2/3 + l)2\ 
2m 
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Similarly, since only a quarter of dj_2 might be in the terrain, the probability of this area is 
bounded from below by 7rrj^_2/(4/3)- Plugging this into the analysis of Lemma [3.1^ we have 



As before, we thus have 



vrr. 



i-2 



4/32 



m— 1 



E 



L j i>2 



< exp 



^ E E 1^(2/3 + 1)" ■ exp(-(2/3 + 1)^-^) = 0(1). 

j=l i>2 



Tir, 



i-2 , 



(m-1) ) < exp(-(2/3 + 1)2^"^) 



The only missing component is bounding the expected number of points of P fl 62, as they 
can affect the Voronoi diagram in □. Arguing as above, this quantity is bounded by m/3 • 
(area of 62) < mPnrl < f (2/3 + 1)^ = 0{I3^). 

Note that if we drop the factors of /3 in the bounds for X/, , and the area of d2, then 
the above becomes a proof for the case when we sample from the unit square ■ 

Theorem 3.3 Let 7 he a terrain. Let P he a random sample of m points, either sampled 
uniformly from the surface of 7, or uniformly from the domain and then lifted vertically up 
to the surface. The expected comhinatorial complexity ofVor{P) is in 0{^/3^X{n + m)) . 



Proof: As described in Section 12731 the combinatorial complexity of the Voronoi diagram is 
the sum of the number of breakpoints of Voronoi edges, the number of Voronoi vertices and 
the number of intersections of triangulation edges with chords of Voronoi edges. By Fact 12.41 
the number of breakpoints is bounded by 0{n), and by Lemma [2.61 the number of Voronoi 
vertices is bounded by 0{m). 

It remains to bound the number of intersections of the set of chords with the triangulation. 
To this end, we place a grid on the domain of the terrain, such that the side length of each 
grid cell is / = I/a/tti and we obtain M = 0{m) grid cells which together cover the domain 
of the terrain. Now, let Ci, . . . , Cm denote these grid cells. Consider the grid cell Cj and 
the set of chords of the Voronoi diagram which intersect this grid cell, let this set be Bj. 
Similarly, let denote the subset of edges of the triangulation, which intersect Cj. Since we 
assumed that the triangulation is A- low density, also Ej is a A-low density set. By Fact 12.51 
we have that the set of chords, which originate from the same Voronoi edge (and therefore 
from the same bisector) form an 0(.^)-low density set. Let ki denote the number of Voronoi 
edges that contribute chords to Bj, we have that B,j is a 0(^/cj)-low density set. By Fact 12.11 
the number of intersections between objects of Ej and objects of Bj is in 0(^/cj |Ej| + A |Bi|). 

Now, in order to bound the overall number of intersections, let Bf ' denote the subset of 
chords which are longer than /, similarly, let B^-' denote the chords in Bj which have length 
smaller or equal to / and let E^' and Ef^ be defined analogously. By the above analysis, 
we have that there exists some constant ci, such that it holds for the overall number of 
intersections Xi 



M M 

X<ci Y^i^h |E.| + A |B,|) = ci ^(eA:.(|E>'| + |Ef |) + a( 

i>l i>l 



B 



B 



<i 
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By the definition of low density sets, we liave tliat |E,f | = 0(A) and |Bf | = O(^fcj), since 
tliey intersect tlie bounding ball of the grid cell Cj, which has radius 0{l). Therefore, it must 
be that there exists a constant C2 such that. 



M 



AI 



+ A 



B 



<i 



C2Hki +ciA4|B| 



where the last inequality follows from the fact that any chord in Bp can intersect at most 
four grid cells, since the grid cells have side length equal to / (similarly any edge in Ep can 
intersect at most three grid cells). 

Finally, note that the number of Voronoi cells that are expected to intersect a grid cell 
in their projection is bounded by 0(/3^) by Lemma [3.21 Thus by Corollary 12 .71 we have that 
E[/ci] = 0{(3^). Therefore in expectation. 




ciA4|B| 



M 



-<l 



+ C2A^E[fci]) +ciA4|B| 



4|B| <C3e/3'A(3|E| 



M + 4|B|) 

for some constant C3 (note that we used the fact that Ep is independent of the random 

sampling). Furthermore, observe that by Lemma [2.61 the overall number of Voronoi edges 
is 0{m). Recall that every Voronoi edge is broken up by breakpoints into chords. Every 
breakpoint increases the number of chords by one. Using Fact 12. 4[ it follows that the overall 
number of chords |B| is in 0(?t, + m). Since the number of edges of the triangulation |E| is 
0{n), we conclude that E[x] = 0(^/3^A(n + m)). ■ 



4 Lower bound 

In this section we show that if we drop the assumptions on the terrain, then the expected 
worst-case complexity of the resulting geodesic Voronoi diagram can be f2(nm^/'^) if the sites 
are sampled uniformly at random from the unit square. 

In the following we will refer to the walls of the unit square defined by x = 0, x = 1, 
y = 0, and y = 1 as the west, east, south, and north walls, respectively. 

Example 4.1 We start with the easy case of planar map and points in the unit square. 

First place m points in a column near and parallel to the west wall 

of the unit square such that the spacing between each adjacent pair of 

points is 6(l/m) (we assume m = 0{n)). The planar map consists of n 

vertical lines near the east wall of the unit square that extend from the 
north wall to the south wall. Now, the boundaries of the Voronoi cells 

of these points extend from the west to the east wall and are parallel 

to the north and south walls, and hence the complexity of the overlay ' LLLLLLLL 

of the Voronoi diagram with the planar map is 0(nm), since it is an n x m grid. 

The later constructions use this as their starting point. 
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4.1 Farming — an Q{n^/rn) example 



4.1.1 Construction 

The height function used in the following construction of a terrain has (essentially) only two 
values, zero and h. The areas between a part of the terrain that is of height zero and of 
height h consist of very narrow and steep boundary regions. This intermediate boundary 
would have a very small measure in the projection, and the reader can think of it as having 
measure zero. Moreover, h is chosen to be sufficiently large so that no point at height zero 
can affect the Voronoi diagram at height h. One can therefore view the following terrain 
construction as a flat unit square, where we have cut out or "forbidden" areas (that have 
height 0). Therefore, for the sake of simplicity of exposition, an area being constructed is 
fiat, at height h, and the adjacent forbidden area is at height zero. Our main building blocks 
will be farms. We define a farm to be a square of side length l/{c^/rn). Intuitively, farms 
are part of the terrain which with constant probability (the constant will depend on c) will 
receive at least one point from the random sample (i.e. a farm takes the place of a point 
from the example in Exercise 14. ip . We define the diameter of a farm to be the quantity 
5 = ^J^jic^fm) (that is, the distance of the furthest two points in a farm). 

We now define a sequence of ridges to take the place of the planar 
map from Exercise 14.11 (i.e. in expectation we would like the Voronoi 
diagram to look like a grid over the ridges). Formally, let a sequence 
of ridges of length be a sequence of 2n rectangles, ri, . . . , such that the right edge 
of Tj is the same as the west edge of rj+i for i = 1, . . . , 2n — 1, rj has a slope of 45° for odd i 
and —45° for even i, all rectangles extend from the north to south walls of the unit square, 
and the geodesic distance from the left edge of ri to the right edge of r2„ is 1/ {c2"') (which 
is 0(1/2"^) since we assumed m = 0{n)). Refer to the figure above to the right for a side 
view of the ridges. 

The construction of the ^}{n^/m) example is as follows. Place 
G(i/rn) farms from north to south along the west wall of the unit 
square, with 2/(cA/m) spacing in between each adjacent pair. Next 
build a sequence of 6(n) ridges near (and parallel to) the east wall of 
the unit square. Then connect each farm directly to the leftmost ridge 
by creating a line parallel to the north and south walls connecting the 



AAAAA 



Side view of ridges 



south-east corner of the farm to the first ridge. See figure on the right. We refer to such a 
line as road. The roads stay at height h and to the left and right of a road, the height drops 
to zero as described earlier. 



4.1.2 Analysis 

Definition 4.2 The point at which a farm connects to its road is its entrance (i.e. the 

south-east corner of the farm), and the point at which the road connects to the leftmost ridge 
is its exit. We say that the point (from the random sample ofm points) that is closest to the 
entrance for some farm, is that farm's dominating point. Let p and e be the dominating 
point and exit, respectively, of some farm. We say that another point q from the random 
sample that is contained in another farm and such that djlgjc) < d'j(p, e), eliminates (the 
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Voronoi cell of) p, where dj{p,e) denotes the shortest path on the terrain from p to e. If 
there are no points which eliminate a given dominating point, then the dominating point is 
alive. 

Lemma 4.3 For the construction of the terrain described above, if one picks uniform at 
random m points in the unit square, their induced geodesic Voronoi diagram on this terrain 
has complexity ^}{n^/m). 

Proof: The area of each farm is G(l/m) and hence a sample of m points picked uniformly at 
random from the unit square, will have at least one point with constant probability in each 
farm. Moreover, since we constructed Q{\/m) farms, this implies that in expectation Q{y/rn) 
farms will receive at least one point. Furthermore, the width of the sequence of ridges and 
roads was chosen such that the probability that either receives a point is exponentially small 
(and hence in the following we assume they do not receive any point). 

Now consider a farm which received at least one point, and let p be its dominating point. 
Observe that the Voronoi cell of p contains the entire road connecting this farm to the ridges, 
and its Voronoi cell extends all the way to the rightmost edge of the sequence of ridges, and 
hence will be of complexity Q{n). Indeed, by our construction, only a point from another 
farm can prevent the Voronoi cell of p from reaching the rightmost ridge. However, the 
spacing of the farms was chosen to prevent this. In the worst case p is in the north-west 
corner of its farm, and an adjacent farm has a point q at the south-east corner. Let / be 
the length of a road. Let Cp (resp. Cq) be the exit of the farm containing p (resp q). Now 
consider the geodesic shortest path connecting Cp to the rightmost ridge. Every point on 
this segment is in distance at most 6 + 1 + l/(c2") from p. However, the closest point on this 
segment from g is at a distance of at least / + 2/{cy/m) > 6 + I + l/(c2"), and so q cannot 
prevent the Voronoi cell of p from reaching all the way to the rightmost ridge. 

Therefore, in expectation, we have that Q{^/rn) farms have a point whose Voronoi cell 
extends all the way across the sequence of ridge, which gives a Voronoi diagram that in 
expectation has complexity f^,{n^/rn). m 

4.2 Industrial farming — an Q(jim'^^^) example 

The challenge in improving the example above is that the distance of a dominating point 
to the exit of a farm has too much variance (i.e., ^1/m). Since there does not seem to be 
a way to decrease the variance directly, instead we connect all the farms to the ridges, and 
carefully argue about the expected complexity of the generated Voronoi diagram. 

4.2.1 Construction 

In the following we assume that m = 0{n). 

We set the side length of each square farm to be Xj ^fm and construct a sequence of Qin) 
ridges near the east wall of the unit square. We will place an M x M grid of farms inside the 
unit square, where M = \_^fmj^\. Specifically, the spacing between columns (which extend 
from north to south) will be 1/ ^Jm and the spacing between rows (which extend from west 
to east) will be Ij^pm. The grid starts in the north-west corner of the unit square. 
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We now describe the connecting roads from the farms to the ridges. The following 
construction of the roads will ensure that the length of each road is the same and that the 
distance between adjacent exits on the ridges is at least 1/m. These two properties will be 
sufficient for the analysis in the next section to go through. 

Consider a given row of farms. Number the farms in this 
row /o, . . . , /m-1 in increasing order of their distance to the 
west wall. Every farm has dimensions x and 

the spacing between two consecutive farms in a row is 1/ y/m. 
As such, the x coordinate of the entrance of the ith farm is 
Xi = {2i + 1)/ \/rn (as before, the entrance to each road will be 
at the south-east corner of the farm). The directions the ith 
farm's road goes from entrance to exit is described as follows: 

(a) south for a distance of at^a = i/m, 

(b) west for a distance of = {xi + ai^a)!'^, 

(c) south for a distance of ai^c = l/y/rn — 2ai^a, and 

(d) east for a distance of ai^d = w — (x, — ai^h) all the way to the first ridge, where w is 
the distance from the west wall to the first ridge. 

This layout is sketched in the figure above to the right. Note that the spacing in this 
figure only approximately matches the description. 

Sanity checks. The road of the zth farm starts at x coordinate Xj, goes west for a distance 
of aj^b and east for a distance of ai^d- Observe that the x coordinate of the exit of this road 
is Xi - ai^b + oti,d = Xi- Ui^b + w - {xi- Ui^b) = w. 

Observe that the ith farm will connect to the ridges in north to south distance aj,a + aj,c = 
1/ y/m — i/m from the southern boundary of the row of farms. That is, adjacent farms in the 
row have exits in distance i/m apart along the first ridge (exits are y/m) apart between 
rows). Furthermore, each road is of the same length. Indeed, let be the length of the road 
for the zth farm to its exit. We have that 

, Xi -\- Q-i a ,1 „ , ( Xi -\- Q-i a \ 

Ti = ai^a H ^ h ^= - 2ai^a + w -IXi 

2 y/m \ 2 y 

^ V ' " V 

That is, all the roads have exactly the same length. 

4.2.2 Competing farms 

We now prove that in expectation 0(m^/^) dominating points will have their Voronoi cells 
reach all the way to the east wall across the sequence of ridges. 

Observation 4.4 Let p and Cp be the dominating point and exit, respectively, of some farm 
/. Let g be a point which is in some other farm /' with exit Cg. If /' is i farms away (in 
the north to south order of the exits of the farms along the first ridge) then ||ep — eg|| = i/m 
and hence q eliminates p if and only if (^^{q, Cp) = dT(g, Cg) + i/m < dj{p, Cp). 
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Next, we prove that each dominating point is ahve with probabihty f2(l/m^/^), using 
the following helper lemmas. 

Lemma 4.5 Let X be a positive random variable with expected value /i. We then have that 
E[e-^] > 6-2^^/2 

Proof: Markov's inequality implies that Pr[X < 2/i] = 1 - Pr[X > 2/i] > 1 - fi/2fi = 1/2. 
Therefore, by the definition of expectation, we get E [e-^] > Pr[X > 2^] * + Pr [X < 2/i] * 
g-2/. > e-2^/2. ■ 

Lemma 4.6 Let f be a farm, and let r{f) be the random variable that is the distance of the 
closest site (that falls into this farm) from the farm entrance (if there is no site in this farm, 
we set r{f) = ^J^jm). Then, for any distance s, we have Pr[r(/) < s] < ms'^ir/A, where 
equality holds for < s < 1/ ^/m. 

Proof: Recall that the entrance of a farm is at the south-east corner. Hence a point in the 
farm which is in distance at most s from the entrance must fall into the intersection of the 
farm with a circle of radius s whose center is at the entrance of the farm, see figure. 

Therefore, if the radius of the circle is less than the side length of the 
square, i.e. s < iZ-y/m, then the intersection is a quarter disk and so the 
area is exactly tts'^/A. Otherwise, the top and left portions of the quarter 
disk will be needed to be clipped to the farm and so the area is < irs'^/'i. 
Now, the probability of the ith site to fall into this disk is < irs'^/'i, and 
since we sample m sites (independently), by the union bound the claim 
follows. ■ 

Lemma 4.7 Let p and e be the dominating point and exit, respectively, of a farm f . Let 
r = da-(p, e). Let fi be a farm which is i farms away from f (either in north or in south 
direction), and let Xi be the number of points which fell into fi. Let denote the probability 
that no point from fi eliminates p (see Definition \4.2\ ). Then aXi > exp(— m(r — i/rnf'Xi/2). 

Proof: Let gi, . . . , qxi be the Xi points that fall into farm fi. Let be the exit of /«, and let 
dj = dr{qj, Cj), for j = 1, . . . , Xi. Arguing as in Lemma we have that Pr[dj < s] < s'^ir/A 
for all j. By Observation 14.41 a point qj eliminates p if and only ii dj < r — i/m. For j ^ I, 
whether or not qj or qi kill p are independent events and hence 




oix = I I Pr Ml > r — i/m >(l , ^ 1 > expl —m(r — i/m) 

' fJi ^ J V area(/) J V 

Lemma 4.8 Let p and e be the dominating point and exit, respectively, of some farm f. Let 
r = d-j{p,e). Let Xi (resp. Yi), for i = 1, . . . , \ rm\, denote the number of points which fall 
into the farm which is i farms to the north (resp. south), from f in the order of the exits 
along the first ridge. Then the probability that p is not eliminated given the values of Xi and 
Yi for all i, is at least 

-m ^(r -i/mf7i{Xi + Yi)/2 
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Proof: First note that for i' > [rmJ, we have that i'/m > ([rmJ + l)/m > r. Namely no 
point from a farm i' farms away can kill p, and hence we can ignore such farms. 

Given the value Xj and Yi for all i, whether a farm contains a point which eliminates p 
is independent from whether any other farm contains a point which eliminates p. Therefore, 
by Lemma 14. 7[ 




Xi , . . . , X^rm\ 1 Yi , . . . , ^[rmj 



[rmJ \ / [rm.J 

n n 



[rmJ 



— m r — I 



[rmJ 



/m)VXi/2) JJ exp(-m(r - z/m)27rF,/2) 



exp -m ^(r - z/m)V(Xi + Fi)/2 



i=l 



Lemma 4.9 Let p and e be the dominating point and exit, respectively, of some farm f , and 
let r = da-(p, e). Then the probability that p is alive is > | exp(— 2r'^m^). 

Proof: Let Xi and Yi, for i = 1, . . . , \rm\ , be random variables equal to the number of points 
which fall into the farm which is i farms to the north or south, respectively, from / in the 
order of the exits along the first ridge (note that if there is no farm i farms to the north or 
south, then Xj = or Fj = 0, respectively^ . 



Lemma 14.81 tells us that Pr 

\ [rmJ 



p is alive 



Xi , . . . , X|^j.^j , Yi , . . . , F^rmJ 



> e ^, where T 



Yli=i "~ 'i/m)^7r(Xj + Yi)/2. Since the area of each farm is 1/m, We know that E[Xi] < 1 
and E[Kj] < 1 for i = 1, . . . , [rmJ ("< 1" is used instead of "= 1" since there might not be 
a farm at that distance). Therefore, 



E[T] = E 



[rmJ 

E 

2=1 



m{r - i/mY7r{Xi + Yi) /2 



[rmJ 



^ m(r - t/my7r{E[Xi] + E[Yi\)/2 



i=l 



[rmJ 



1=1 

By Lemma [4.51 Pr p is alive > E[e~^] > -exp(— 2E[T]) > - exp(— 2r^m^) . 
Observation 4.10 Lemma [4.91 implies that if r < then 

Pr[p is alive] > -exp(— 2r^m^) > — -. 



Lemma 4.11 The probability that a farm f gives rise to a Voronoi cell that is not eliminated 
IS > Ti/{%e^nv'/'^). 
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Proof: By the equality in Lemma 14.61 for short distances, we know that Pr[r(/) < s] = 
ms^7r/4, for < s < Xj^Jm. Let p be the dominating point of /. By Observation 14. lOl and 
Lemma 14.61 we have 



Pr[p is ahve] > Pr 



> — Pr 



(r(/) < m ^/^) n(p is ahve) 
r(/) < m-2/3l < ^ 



= Pr 

nm 



p is ahve 

7T 



r{f) < Pr[r(/) < 



2e2 4m^/3 Se^m 



1/3 



Theorem 4.12 /n expectation, the Voronoi diagram will be of complexity f2(nm^/^). 

Proof: Every farm which receives a point from the random sample has a dominating point. 
Since 6(m) farms were built, and each farm receives one point in expectation, the expected 
number of dominating points is 6(m). Therefore, by Lemma [4.1H the expected number of 
alive dominating points is f2(m * (7r/(8e^m^/^))) = r2(m^/^). 

Given that a dominating point is alive, the probability that its Voronoi cell does not 
reach the rightmost ridge is negligible. Therefore, in expectation, if there are f2(m^/^) alive 
dominating points and 0(n) ridges then the complexity of the Voronoi diagram will be 



5 Conclusions 

We investigated the expected combinatorial complexity of geodesic Voronoi diagrams on 
polyhedral terrains in two settings where the sites are being picked randomly. Usually, such 
random settings are the great simplifier - for example, the expected complexity of the convex 
hull of n points picked uniformly in the unit square is O(logn) - but in our case the situation 
is considerably more subtle. 

We proved that the expected complexity is linear if one assumes low density and bounded 
slope and the domain of the terrain is a unit square. On the other hand, we described a 
worst-case construction of a terrain which implies a super-linear lower bound on the expected 
complexity if these assumptions are dropped. This implies that the probabilistic analysis 
alone does not yield a linear complexity. 

There are still many interesting open question for further research. 
In particular, if we relax the realistic input assumptions, is the expected 
complexity still linear or can one show other lower bounds? One could, 
for instance allow the terrain to have a constant number of triangles, 
where the slope is unbounded, or drop the steepness assumption com- 
pletely. Consider the farming layout to the right, where the sequence of 
ridges have been replaced by a recursive low density construction. One 
problem with this construction is that the low density assumption requires the ridges to have 
a non- negligible area, which could catch points from the random sample. 
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A How fat are Random Voronoi Cells? 

In light of the known results on Poisson Voronoi Diagrams lOBSCOOj . the results in this 
appendix are not too surprising. We include our analysis here because it is relatively simple 
and it is self contained. We are unaware of any work analyzing the fatness of Voronoi cells. 
Furthermore, our sampling model is different than the one used in the Poisson Voronoi 
Diagrams. 

Let P be a set of m points picked uniformly from the unit square [0, 1]^. To simplify the 
presentation, we assume the square has the torus topology (i.e., identifying together facing 
edges). Let pi be the first sampled point, and let P = {pi, . . . , Pm}- In the following, let 
cell(pi, P) denote the Voronoi cell of pi in the Voronoi diagram of P. 

Lemma A.l For any integer j > 0, with probability > 1 — Qe^~^ , the diameter of C = 
cell(pi, P) is bounded by Rj = 4 ■ jj {m — 1). 

Proof: Partition the plane around pi into 6 equally spaced cones, of 60 degrees each. Consider 
such a cone, and break it into tiles of area l/(m — 1) each, by cutting it by parallel lines 
forming equal angles with both sides of the cone. In particular, let Tj be the ^th tile in this 
partition, with the tiles ordered in increasing order way from the center pi. 
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The first j tiles of tfiis cone form an equilateral triangle of 
area j/ (m — 1). As such, letting ij and hj be the edge length 
and height of this triangle, we have that hj = {\/3/2)£j, and 
its area is hjij/2 = i'jy/^/A = j/{m — 1), which implies that 

ij = 2- 3-i/Vj7(m-l). 

Clearly, the probability that the first j tiles would not 
contain any point of P is (1 — j/ (m — < exp(— j(m — 

l)/(m — 1)) < exp(l — j). As such, the probability that the 
first j tiles in any of these six cones is empty is bounded by 
6exp(l — j). Namely, with probability > 1 — 6exp(l — j), we have that all cones have a 
point in them in distance at most ij from pi. This implies that the Voronoi cell of Pj in such 
a cone in contained the in first j tiles. As such, the diameter of this cell is bounded by 2ij, 
as claimed. ■ 




Lemma A. 2 Let X be the distance from pi to the second closest point to it in P \ {pi}. Let 
i{m~i)n ■ '^^ indicator variable that is one if and only if X & [rj+i,rj]. For 



i>A, we have Pr[Xi = 1] < 1/r. 

Proof: In the following, let pi = nrf = l/{i{m — 1)) be the probability of a random point 
to fall inside the disk of radius centered at p. The probability that this disk has exactly 
k points in it is exactly 



These probabilities fall quickly with k, and can be treated as a geometric series. Indeed, 
Ofc+i k\{m — 1 — k)\pi (m — 1 



k)pi ^ 2{m - 1 



k) ^1 



ak {k + l)\{m - 1 - k - - Pi) {k + - Pi) - i{m - l){k + 1) - 2' 
for i > 4. Now, we have 



Pr 



X e [ri, ri+i]l < Pr [|disk(p, r^) n (P \ {pi}) | > 2] = ^ < 2^2 



as desired. 



k>2 



< 2 



(m — l)(m — 2) 



1 

< -r. 



z^(m — 1)^ ' 



Lemma A. 3 If Xi = 1 then there is a disk of radius rj+i/4 contained inside cell(pi, P). 



Proof: Let q and s be the closest and second closest nearest neighbor to pi, respectively, in 
P\{pi}- 
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Consider the line £ through pi and q, and place a point t on this 
line in distance rj+i/4 from pi, on the side away from q. We claim 
that the disk D of radius rj+i/4 centered at t is fully contained in 
cell(pi, P). Indeed, observe that the Voronoi cells of all the points 
of P in distance at least rj+i from pi can not contain any point of 
D. As such, the only point that might be problematic is q. But 
then, observe that the bisector between pi and q is orthogonal to 
the line and it can not intersect thus implying the claim. ■ 

The fatness of shape X is fat(X) = R{X)/r{x), where R{x) (resp. r(x)) is the radius 
of the smallest (resp. largest) disk enclosing (resp. enclosed in) X. 

Theorem A. 4 The fatness of C = cell(pi, P) is constant, in expectation. 

Proof: Let Yj be an indicator variable that is one if and only if the diameter of cell(pi, P) is 




in the range 



Rj-i, Rj 



By Lemma [A.lt we have that Pr[^- = 1] < 6e^ ^ ^. Similarly, let 



Xi be as in Lemma IA.2I Note, that if Xj = 1 and Yj = 1 then 
F = fat(C) < 



ri+i/4: v^4/((z + l)(m- l)7r) 



As such, we have E[fat(C)] < J2Zi EJli P^i^i = 1) n {Yj = 1)] 4:^j{i + l). Setting L( 
10 [In z] + 1, we have that 

oo L{i) oo oo 

E[fat(C)]<5^^Pr[X, = l]4V7(7Tl) + 5^ Pr[Y, = 1] 4,/JilTl) 

i=l j=l i=l j=L(i)+l 

oo _ oo oo 

<E^4L(.)v/L(z)(z + l) + X; E 6e^-^-i-4V7(^ 



i=l 
oo 



i=l 



o 



1=1 



i=l j=L{i)+l 



0(1] 
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